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ABSTRACT 



Randomly oriented relativistic emitters in a relativistically expanding shell provides an al- 
ternative to internal shocks as a mechanism for producing GRBs' variable light curves with 
efficient conversion of energy to radiation. In this model the relativistic outflow is broken into 
small emitters moving relativistically in the outflow's rest frame. Variability arises because an 
observer sees an emitter only when its velocity points towards him so that only a small frac- 
tion of the emitters are seen by a given observer. Models with significant relativistic random 
motions require converting and maintaining a large fraction of the overall energy into these 
motions. While it is not clear how this is achieved, we explore here, using two toy models, the 
constraints on parameters required to produce light curves comparable to the observations. We 
find that a tight relation between the size of the emitters and the bulk and random Lorentz fac- 
tors is needed and that the random Lorentz factor determines the variability. While both models 
successfully produce the observed variability there are several inconsistencies with other prop- 
erties of the light curves. Most of which, but not all, might be resolved if the central engine 
is active for a long time producing a number of shells, resembling to some extent the internal 
shocks model. 

Subject headings: gamma-rays: bursts; turbulence; 



1. Introduction 



GRB's temporal variability played a major role in the understanding how GRBs operate. Standard ex- 
ternal shocks, in which the external medium slows the relativistic ejecta, cannot produce efficiently variable 
light curves ( |Sari & Piran||1997[ ). While internal shocks resolve the variability and agree with other prop- 
erties of GRB light curves (e.g. Nakar & Piran |2002b| Ramirez-Ruiz & Fenimore|2000 ) they suffers from 
several drawbacks. First and foremost is their low efficiency (Kobayash i et al.|1997 Daigne & Mochkovitch 
1998] see however, |Kobayashi et al.||1997[|Kobayashi& Sari|2001[|Beloboro dov 2000). This is particularly 
troublesome in view of the high efficiency implied from comparison of the prompt 7-rays luminosity and 
the kinetic energy that remains in the outflow. Detailed models for the emission mechanisms of the prompt 
7-rays pose other problems (Kumar & McMaho n |2008 ). 
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Extemal shocks can produce highly variable light curves if the outflow is slowed down by small external 



clumps. Each clump producing a short pulse. However, this process will inevitably be inefficient (Sari & 



Piran|1997] ) as the overall covering factor of the emitting regions is St/T (5t and T are the pulses' and the 



burst's durations). Observed values of 6t/T are typically ~ 0.01 and in can be as low as 10 4 (Nakar & 
Piran |2(X)2a) . 



Lyutikov & Blandford (2002, 2003) (see also |Lazar|[20031 [L05], [Lyutikov|[2006l |Narayan & Kumar 



2008) proposed that variability can be recovered while maintaining high efficiently it the shell that moves 



with a bulk Lorentz factor F contains emitting clumps (see fig. [T} that move with random macroscopic 
relativistic velocities (with a Lorentz factor 7'). A clump is observed only when its radiation cone (with 
an opening angle of the order of 1 /F^' in the lab frame) points towards the observer. The filling factor of 
the clumps may be unity, recovering high efficiency. However, as only a small fraction of the clumps are 
observed at any given time, the light curve can show rapid variability. The overall duration of the burst, is 
the larger between the angular time and the shell's light crossing time (max{ R/cF 2 , A/c}), where R and A 
are the shell's radius and width respectively and c is the light speed, allowing for emission radii much larger 
than 5tT 2 , required in the internal shocks model. The temporal variability is then dictated by the random 
Lorentz factor, 7', reflecting the activity of the emitting region and not those of the inner engine. 

While it is unclear how macroscopic random relativistic motion can be generated, we assume that it 
does and examine, using two simplified toy models that includes the essential ingredients, the conditions 
under which the temporal features of the observed light curve can be produced (see L05). We describe our 
first toy model, which we call here relativistic turbulence, and derive analytic constraints and numerical light 
curves in £j2| In ^3] we consider a second toy model proposed by Lyutikov (2006) that is based on sub-jets 
and compare it with the first one. We summarize the results and compare both models with observations in 





2. Relativistic Turbulence 



Our (L05) kinematic toy model for relativistic turbulence considers a shell which is divided into discrete 
randomly distributed emitters that have randomly oriented relativistic velocities. The emitters change their 
direction of motion continuously, as expected in a turbulent medium. Since the emitters exhibit a coherent 
macroscopic motion, we require that each emitter is causally connected, and that it changes its direction 
on times longer than the causal time scale but shorter than the shell crossing time. The length scale of the 
emitter in its own rest frame, I", is assumed to be similar in all three dimensions, and it emits isotropically 
in this frame. We define a dimensionless scale tp = I" /R, which is the angular scale of an emitter that 
points towards the observer (see inset of fig[T]). The emitters radiate as the shell moves from Rq to 2Rq. 
Due to the turbulent motion the positions and directions of the emitters change with time. We model this 
by a set of successive shells between Rq and 2Rq. Each new shell is constructed with randomly distributed 
emitters, representing the random changes in direction of the turbulent motion. The time difference between 
two shells is, r', the time it takes for the emitters to turn an angle of 7 /_1 (in the shell frame). 
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Turbulence Sub-Jets 



Fig. 1. — The basic kinematic model - relativistically expanding shell with an ensemble of emitters, which 
are moving at random relativistic velocities within the shell. The inserts describe the geometrical details 
and definitions of the two alternative models. Left: relativistic turbulence - emitters of size Rip move 
relativistically; Right: sub-jets - energy from regions of size Ripsj is extracted into relativistic jets. 

Note that there are three frames: The lab frame; the shell's frame, denoted by a prime, which is boosted 
radially with a Lorentz factor T relative to the lab; and the frame of each emitter, denoted by two primes, 
which is boosted by (randomly oriented) 7' relative to the shell frame. The observer is, of course, at rest 
relative to the lab frame. However, the observer time, namely the arrival time of photons (denoted t) differs 
from the lab time by the usual time of flight arguments ( |Rybicki & Lightman |1979 ). 

The Doppler shift from an emitter is: 

A = [7(1-/3 -cos a)}-\ (1) 

where 7, (3 and a are the Lorentz factor, velocity and the angle between the velocity and the line to the 
observer (both in the lab frame). The flux that reaches the observer from this emitter is: 

F u = j J^A^O^J^A 3 ^, (2) 

where /"„ is the specific intensity and the second relation holds for a small enough emitter (D is the distance 
to the observer). An implicit K correction arises from the difference between the v and v" . 

Define 9 as the angle between the line that connects the origin and the emitter and the line that connects 
the origin and the observer (see inset of fig [I]). The maximal Doppler boost, A max = 47T, is obtained 
for emitter at 6 = that moves along the line of sight. The flux decreases like A 3 plus a K correction. 
Therefore, we consider emitters only if A, > A max /2. What is then the probability that an emitter at angle 
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6 is observed, dP/dO = 5(0, r,Y)? At 9 = it is 5(0, r,V) ps l/4 7 /2 while 5(1/T,r,7') = 0. This 
suggests that 5 scales as S(9, T, 7') = y~ 2 S(9T) (L05), implying that the average probability that an 
emitter will be visible from an arbitrary position on the shell is: 

The factor 0.3 was evaluated numerically (L05) and can be ignored at the accuracy level of our discussion. 

The arrival time from an emitter at R, 9 is: 

R-Ro R9^ X 
2cT 2 + 2c c ' W 

where x{< A) is the distance of the emitter from the front of the shell. As the last photons will arrive from 
2Ro, an angle of 1 /T and x = A, the overall duration of the burst will be a function only of A and T (and 
not 7'): 

where we define d = AT 2 /R. As the shell is expected to expand relativistically in its own framed > 1. 
For d > 1 the shell's width, as well as T, are determined by the engine activity time while for d = 1 they 
don't. 

The duration of a pulse arriving from a single emitter is the longest of the three following time scales: 

(i) The duration over which the emitter points towards the observer, namely the duration over which 
the direction of motion varies by an angle 1 /r 7 ' in the lab frame (l/ 7 ' in the shell's frame). As the emitter 
is confined to the shell it should make at least a 7r/2 turn during A'/c, implying that the time to turn by 
l/ 7 ' (shell's frame), r', is shorter than A' /cy'. Causality puts a lower limit on r' of Rtp/c. Therefore: 

R^j/c < t' < A' /cry'. (6) 

In the observer's frame, this translates to: 

Ri/j/T^c < r < A/cV 3 . (7) 



(ii) The emitter's light crossing time in the lab frame (in the direction along the line of sight). For 
an emitter moving towards the observer this time is Rtp/^/T. (iii) The angular time scale - At the largest 
possible angle, where the emitter is still visible by the observer, I/7T, the time difference between the first 
and the last photon would be -iJ^sin (I/7T) w Rifj/cr/T. Overall (ii) and 



(iii) are of the same order and much larger than (i). Thus: 

St « Rtp/cj'T. 



(8) 



1 Note that for a hydrodynamic external shock d < 1 



Sari & Piran 



1997 1 but this might not be relevant here. 
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Using Eqs. [5]and[8]we express, N p , the number of (possibly overlapping) pulses expected in a burst: 

T di 

N p = n p - = n p — (9) 

where n p is the occupation number of pulses at any given observer time (i.e., n p ~3> 1 implies many over- 
lapping pulses while n p <C 1 implies long quiescent periods between isolated pulses). 

The number of emitters is 4irR 2 A'/ (Rip) 3 = At: AT /Rip 3 . The emitters obtain new random directions 
(which differ by more than 1/7', in the shell's frame, than the previous ones) after a time r'. Thus, the total 
number of independent emitters, N to t, is larger by the factor R/(cTt'), the ratio of the total duration over 
which the radius doubles and t'. Finally we introduce a filling factor f < 1 allowing for the possibility 
that not all emitters are active all the time or that space is not fully covered by emitters (/ <C 1 is strongly 
disfavored as the efficiency is always smaller than /). Overall we find: 

4-irf d R 



The condition Np = PN to t yields: 
and using [7] 



fd R do 



7 /3 r 3 ^ 2 ct' ' 



f f 

<n p < -. j± (12) 



d(jTip) 2 ~ p ~ ( 7 Ty>) 3 ' 



We demand n p 1 since many overlapping pulses reduce the observed variability, whereas very 
frequent long quiescent times between pulses are not observed. If the shell is in the freely expanding phase 
(i.e., d 1) n p will be of order unity if: 

7'1 



where k is between 2 and 3. Narayan & Kumar (2008) have pointed out that ip = I/7T if one requires 



that the emitters are of the maximal causally allowed size. Note that n p depends quite sensitively on 7'r^ 
and it increases rapidly if ip is smaller than I/7T. This implies, for example, that a significant number of 
small eddies, which may arise in a turbulent cascade, may be problematic. Using the relations 13 and[9} and 
assuming the causal limit for t'\ 

'^(L) 1/e pL, ( i4) 

\n p J V dSt' 

leading to 7' 10/ \fd for typical values of T/dt. Note that while the model determines 7' it does not 
constrain T and R. 

Fig. [2] depicts simulated light curves (L05) for four choices of parameters. The two upper panels 
have n p = 1 with different emitter sizes. Both light curves are highly variable and densely filled with non- 
overlapping pulses. However, as d = 1, the underlying overall envelope of the pulses is seen. As the emitters 
are smaller on the right panel it has more pulses than the left one. The envelope is observed since only a 
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small fraction of the volume and hence fewer pulses are seen early on. Similarly at t > (d + 1)Rq/2cT 2 
pulses from small 6 values are not seen, implying that only lower amplitude pulses (on average) are observed 
during the last T / (d + 1) of the burst. The envelope is stretched on bottom panels where d = 10. The lower 
left panel depicts a very low n p with a rather sparse light curve. The lower right panel depicts a light curve 
of a wide shell and n p = 0.7, which is rather similar to observed bursts. For n p 3> 1 ( not shown ) the pulses 
are overlapping and all variability is erased, leaving only the envelope. 
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Fig. 2. — Numerical monte Carlo simulations of light curves for different combinations of parameters 
(shown on each frame; the scalings used eliminate the dependence on R and T). The flux of each pulse is 
calculated assuming that the radiation efficiency is constant per unit mass for all emitters in their rest frame, 
namely that /"„ oc (ipR)~ 3 , therefore F u oc A 3 /(ipR). We approximate each pulse as a Gaussian with the 
above parameters and we sum over all contributions to construct a light curve. In both upper frames n p = 1 
and d = 1. In the top left, the emitters are as large as causality allows while on the top right they are 1/2 of 
this value and the filling factor is lowered to compensate. In both the overall envelope is seen clearly. The 
sparsity of pulses is apparent when n p = 0.01 (bottom left) and the "straightening out" of the envelope is 
clearly seen (bottom two panels) when d = 10. 



3. Relativistic Sub- Jets 

Motivated by reconnection in highly magnetized outflow Lyutikov (2006) considers a model in which 
relativistic sub-jets (SJs) are accelerated to a Lorentz factor 7' by dissipation of the bulk energy in many 
different "mini-engines" within the relativistically expanding shell. These "mini-engines" or acceleration 
sites correspond to reconnection sites within the magnetized flow (e.g., |Lyutikov & Blackman 2001 ). The 



"mini-engines" are at rest in the shell frame. Each mini-engine operates for a time t' SJ . The directions of 
the accelerated sub-jets are random in the shell frame but the opening angle and direction of each is constant 
while its mini-engine is active. 
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The sub-jet extracts energy from a region of size I' = 5t' SJ (3c, where f3c is the speed of extraction 
of energy from the surrounding region (relativistic reconnection suggests 0.1, Lyubarsky||2005 ). The 
observed duration is: 

«5j - T ^T' (15) 



where we define the dimensionless parameter ipsj = Slightly generalizing Lyutikov (2006 1 we write 
the probability to observe an emitter as <p 2 /47rT 2 , where cj> = max(\/^sj, 1/tO an( ^ ^SJ * s ^ e surj -jet 
opening solid angle. Following the notation of £j2| we define, /, the filling factor of regions from which 
energy is extracted into the sub-jets (N tot = f47rA'R 2 /l' 3 ). The occupation number of observed pulses is: 

n P*-> = a jt r2 • (16) 

The condition n PSJ rj 1 yields: 

jlvar' (17) 

^WJ- 1 ^- (18) 

This implies that an efficient (/ ~ 1) highly variable burst requires either a large 7' or a wide shell (for 

T/st ~ 100, yd ~ 100/3- 3 / 2 .). 

This constant direction of the emitters and the fact that causality in the shell's frame determines the 
sub-jet size, cSt', are the main kinematic differences between the sub-jet and the turbulence model (in which 
the emitter's direction varies and causality in the emitter's frame determine its size, I"). For the same n p and 
St/T the two models give similar light curves. In particular, an overall (rising and falling) envelope for the 
light curve is expected in the sub-jet model as well. 



and 



4. Discussion and Conclusions 

We have derived conditions on the parameters of relativistic random emitters needed for producing 
variable GRB light curves. This is characterized by n p 1 which ensures that typical pulses don't overlap 
and are not too sparse either. Our numerical simulations show that for 0.03 < n p < 3 one obtains light 
curves that resemble observed GRBs (see fig. [2]). The resulting light curves do not change qualitatively 
when we introduce a distribution of turbulent Lorentz factors and sizes. 

Causality suggests, for relativistic turbulence, that the relation ip = 1 /Tj' between the angular size of 
the emitters, tp, and the turbulent and the bulk Lorentz factors holds naturally (Naraya n & KumaF]|2008 1. 
But, this condition holds when the turbulent eddies are of the maximal possible size and may be broken by 
cascade to lower scales. The condition, ipsj = vT^/m^T arises in the sub-jets model. For high efficiency, 
negligible sub-jet opening angle assuming ~ 1 this reduces to ipsj ~ 1/1^7' or t' SJ R/cTj'. While 
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this is similar to the one obtained in the turbulent model, here there is no apparent physical motivation for 
proportionality between t' SJ and 1/7' and this requires an ad hoc fine tuning. 

In both models the light curves arising from a single expanding shell with d ~ 1 shows a rising and 
falling underlying envelope. Furthermore, a single shell cannot produce bursts which depict long quiescent 
periods. The envelope can be erased if d S> 1, while quiescent periods require an outflow of several shells 
(where naturally d 3> 1). These solutions become marginal in the turbulent model if t' is determined by 
causality, since 7' 3> 1 requires d < 10 (see Eq. 14 1. The sub-jet model, however, may favor d S> 1 as it 
reduces the required value of 7'. 

It seems that with proper conditions (and rather reasonable in the case of the relativistic turbulence) 
these models can produce (efficiently) the observed highly variable GRB light curves. We turn now to 
several shortcomings. First and foremost is the question how such macroscopic relativistic motions can be 
generated and sustained. One needs to convert ~ (1 — 1/7') / of the initial total energy to the kinetic energy 
of the emitters and further dissipation in the emitters' frame is needed to generate the radiation. Additional 
questions involve the shape and other properties of individual pulses versus those seen in observed pulses: 



(i) GRBs show a clear difference between the fast rise and the slow decline of individual pulses (Norris et 



al.|1996 1. The light curve of an individual pulse results from a combination of the motion of the emitter, its 



orientation relative to the observer, its width as well as intrinsic inhomogeneities within the emitter. In the 
relativistic turbulence model the emitter was radiating long before its velocity pointed towards the observer 
and it continues to emit long after it moves away from the observer. There is no reason (on average) for a 
difference between the rising and falling phases of an individual pulsd^ This is not a problem in the sub-jets 
model in which the onset of the pulse corresponds to the beginning of the activity of the emitter. 



(ii) The temporal structure of the first and second halves of GRB light curves are similar ( Ramirez-Ruiz & 



Fenimore j2000| ). The light curves produced in the two models have an overall envelope that favors stronger 



pulses earlier and weaker ones later. This might be resolved by a combination of several emitting shells or 
with very wide shells, but here fine tuning is required in the turbulent model in order to keep 7' 3> 1. 

(iii) Weaker and denser pulses (arriving from emitters not moving directly towards the observer) continues 
at t > T producing the typical envelope of high latitude emission ( [Kumar & Panaitescu 20001. This is 
consistent with rapid declines seen in some the early afterglows. However in many cases the decline is 
faster. In the internal shocks model this is attributed to the dominant contribution of a late pulse, that shifts 
the zero point of the time. Such an option does not arise here unless, once more, we allow for several shells 
or a single wide shell. 

(iv) The duration of an observed pulse is correlated with the preceding interval (Nakar & Piran 2002b[ 
Quill igan et al.|2002] ). There is no reason that such correlation should appear in both models. 



(v) These models predict a Doppler induced correlation between the intensity and E pea k- While stronger 
peaks are typically harder, it is not clear whether this specific relationship is satisfied. 



2 Note that systematic variation of the emitter properties on a time scale of r' will result in a strong signature differentiating 
between early and late phases of the overall light curve, which is not observed. On the other hand non-systematic variations (e.g., 
deceleration and acceleration) are expected to result in similar affects on the temporal structure of rising and decaying parts of 
pulses. 
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We could not find obvious modifications that will address all these issues. While it is not clear that 
those cannot be found, this suggests that the simple versions of these models might not be enough. A simple 
extension of a wide shell d ^> 1 or several separated shells might resolve some of the issues and it might be 
essential for the sub-jet model allowing moderate values of the sub-jet's Lorentz factor. 

This research is supported by the ISF center of excellence in High Energy Astrophysics (TP & AL), 
Marie Curie IRG grant (EN), advanced ERC excellence award and the Schwartzmann chair (TP). 
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